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Abstract. We present a detailed non-spherical modeling of dark mat- 
ter halos on the basis of a careful analysis of state-of-the art N-body 
simulations. The fitting formula presented here form a complete and ac- 
curate description of the triaxial density profiles of halos in Cold Dark 
Matter(CDM) models. This modeling allows us to quantitatively discuss 
implications for shape observations of galactic and cluster halos. The pre- 
dictions of the concordance ACDM model are confronted with the shape 
observations from the Milky Way to X-ray clusters. 



1. Introduction 

It is well known that dark matter halos are not spherical as predicted by current 
theories of structure formation (V. Springel, this volume, for a review). The 
shape of halos depends on the nature of dark matter and cosmological parameters 
in theory, and can be measured through a wide range of observations from the 
kinematics of satellites and debris in galaxies to the hot gas distribution of 
galaxy clusters. In this talk, I will first present an accurate model for the shape 
of dark halos in Cold Dark Matter (CDM) dominated models, and then discuss 
its implications for and confrontations with observed shapes of galactic and 
cluster halos. 



2. Modeling the non-spherical density profiles of dark matter halos 



We use two sets of state-of-the art simulations for the current purpose. The 
first is a set of cosmological N-body simulations with = 512^ particles in a 
100/i~^Mpc box (Jing & Suto 2002), and the other is a set of high-resolution 
halo simulation runs (Jing & Suto 2000). 

To model the shape of dark matter halos, we first find the iso-density sur- 
faces. This begins with a computation for the local density at each particle's 
position. We adopt the smoothing kernel widely employed in the Smoothed 
Particle Hydrodynamics (SPH, Hernquist & Katz 1989). Thirty two nearest 
neighbor particles are used to compute the local density. The left panel of Fig- 
ure 1 shows a typical example of the halo isodensity surfaces. This plot clearly 
suggests that the isodensity surfaces can be approximated as triaxial ellipsoids. 

The isodensity ellipsoids at different radii are well aligned, and the axial 
ratios of the ellipsoids are nearly constant. These facts suggest that the inter- 
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Figure 1. Left panel - A typical example of the projected iso-density 
surfaces within a dark halo. The bottom panels show the triaxial fits to 
five isodensity surfaces projected on those planes. Right Panel - Ratio 
of the quadrupole moments defined in the triaxial model (Qe) and in 
the spherical model (Qs) for five shells at radii from O.OSrvir to O.GSrvir- 
The figure clearly shows that the triaxial model describes the density 
profiles of halos more accurately than the conventional spherical model. 
From Jing & Suto (2002). 
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nal density distribution within a halo can be approximated by a sequence of 
the concentric elhpsoids of a constant axis ratio. To show this to be an im- 
proved description over the conventional spherical description, we compute the 
quadrupolc of the particle distribution within a spherical shell {Qs) or an ellip- 
soidal shell (Qe)- If the spherical (triaxial) model is exact, Qs (Qe) vanishes, 
thus a more accurate halo density profile should yield a smaller quadrupole. The 
ratio of these two quantities is plotted in the right panel of Figure 1. It clearly 
shows that the triaxial model works much better than the conventional spherical 
model. 

Wc found that the matter density within the isodcnsity surfaces changes 
with the major axis R in a way similar to that of the spherical model. The 
density profile can be approximately described by the NFW-like profiles, though 
the concentration parameter is different from the spherical model. In Jing Sz 
Suto (2002), we have given a recipe, based on our analysis of the simulations, 
for predicting the triaxial density profile for a halo in a general CDM cosmogony 
if the halo mass and halo shape (i.e. the axial ratios) are known. 

The mass distribution of CDM halos is well described by the Press-Schechter 
formula (or its modified form) for theories of formation. Thus, in order to 
predict the density profiles for a sample of CDM halos, we need to quantify 
the distribution function p(a/c,b/c) of the axial ratios a/c and b/c at a given 
mass, where a, b, and c are the minor, middle, and major axes respectively. The 
function can be written as 



The probability function p(a/c) and the conditional probability p{a/b\a/c) have 
been measured for halos in our SCDM and ACDM simulations at different 
epochs. An example is given in Figure 2. Jing &; Suto (2002) have found 
accurate universal functions to describe these two probability functions. The 
shapes depend both on the halo mass (through the ratio of the halo mass to the 
non- linear mass M*) and on the cosmological parameters, and the distribution 
of the shapes is very broad. This means that in order to test the shape predic- 
tion of CDM models definitely, one generally needs a large sample of objects, 
e.g. galaxies or clusters of galaxies. As a reference, it would be interesting to 
have a look at typical values of halo shapes. For the concordance ACDM model 
of a fluctuation amplitude = 0.9, the typical values of the axial ratios a/c 
and b/c are 0.7 and 0.8 for Milky Way halos and 0.45 and 0.65 for clusters of 
galaxies. The galactic halos at redshift z = 0.5 are expected to have a/c = 0.5 
and b/c = 0.7, much more eccentric than their counterpart at 2; = 0. 

3. Comparison with the observations 

From the kinematics within galaxies or from the X-ray observations of clusters, 
one may measure the potential of their dark halos. For the triaxial model in 
the last section, Lee &; Suto (2003) presented a calculation for the distribution 
of iso-potential surfaces. Using the eccentricity e = (1 — a^/c^)-*^/^, the ratio of 
the halo potential eccentricity to the halo density eccentricity at radius r is a 
function of r/Ro (Lee &: Suto 2003; the left panel of Figure 3), where Rq is a 



p(a/c, b/c)d{a/c)d{b/c) 



= p{a/c)d(a/c) p(b/c\a/c)d{b/c) 
= p{a/c)d{a/c) p(a/b\a/c)d{a/b) . 



(1) 
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Figure 2. Left Figure- The distribution of the axis ratio a/c of the 

halos in the cosmological simulations of the ACDM model before {left 
panels) and after {right panels) the scaling is applied (see Jing &; Suto, 
2002). Solid, dotted and dashed histograms indicate the results for 
halos that have M4 x 10^ particles within the virial radius. Right 
Figure - The conditional distribution of the axis ratio a/6 of the halos 
in the cosmological simulations of the ACDM model for a given range 
of a/c. The smooth solid curves are the fitting formulas. Prom Jing & 
Suto (2002). 



scale radius of the major axis in the generalized NFW form (Jing & Suto 2002). 
The potential becomes rounder at an outer radius, because the halo mass is 
concentrated in the central region. For a galactic halo of Milky Way mass, the 
axial ratio a/c of the potential is 0.85 at the central region and is larger than 0.9 
at the virial radius. These values are well consistent with the potential shape 
measured for the Milky Way (Sackett, P., Sellwood, J., this volume). 

X-ray morphology of clusters can be computed from the potential under the 
hydrostatic equilibrium assumption, or can be measured directly from hydro/N- 
body simulations. Flores et al. (in preparation) have carried out a comparison 
of the cluster ellipticity between the concordance ACDM model and an Einstein 
X-ray cluster sample. The distribution of the model clusters is completely consis- 
tent with that of the observed clusters, though larger samples both of simulation 
clusters and of observed clusters are badly needed for an accurate assessment. 
Flores et al. also noted that the shape distribution found in their work is well 
consistent with Jing &; Suto (2002). In a recent work of Plionis (2002), a signif- 
icant evolution is found of the cluster X-ray shapes, as higher-redshift clusters 
are more elongated than the local clusters. Floor et al. (2003) found it diffi- 
cult to explain this ellipticity evolution within the concordance ACDM model. 
I would like to point out, however, the comparison between the model and the 
observation carried out by Floor et al. might not be proper, because they used 
individual cluster simulations in their study. The "clusters" they identified at 
high redshift must be smaller than the clusters at z = 0, while in the observation 
of Plionis (2002) clusters at high redshift are much more massive than those at 
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Projected shape of galactic halo 




Figure 3. Left Figure- The ratio of the eccentricity of the iso- 
potential surface to that of the halo iso-density surface from the per- 
turbative result for the NFW profile (a = 1) and for the Moore et al. 
like profile (a = 3/2, dashed line). From Lee & Suto (2003). Right 
Figure - The ellipticity distribution of projected galactic halos at z = 
and z = 1 in the concordance ACDM model. 



z = 0. With the triaxial model of Jing & Suto (2002), it is not difficult to ex- 
plain the redshift dependence of the cluster ellipticity found by Plionis (2002), 
at least qualitatively, because the high redshift clusters in the observation are 
more massive. 

It is promising to probe the halo shapes directly with the weak lensing effect. 
The weak lensing measures the projected mass distribution (Schneider, P., this 
volume; Gavazzi et al. 2003). We have measured the shapes for the projected 
halo mass distribution in the simulations. The right panel of Figure 3 shows the 
distribution of the projected ellipticity for galactic halos (mass about lO^^M©) 
in our simulations. The halos have a mean ellipticity e = 0.31 or axial ratio 
a/b = 0.72 at redshift z = 0, and e = 0.42 or a/b = 0.64 at z = 1. Recently, 
Hoekstra et al. (2003) have tentatively measured the shape of galactic halos 
through weak galaxy-galaxy lensing, and found a/b = 0.66 for galactic halos at 
z ~ 0.5. Their results are in good agreement with our theoretical predictions, 
though it is important to note that they used a weighting procedure to measure 
the halo shapes that must be properly modeled in future comparisons between 
the observation and the theoretical predictions. 

In addition, we have compared the shape distribution of rich clusters in the 
Standard Cold Dark Matter (SCDM; Qq = 1 and Aq = 0) and in the concordance 
CDM (r^o = 0.3 and Aq = 0.7), and found their difference is not dramatic. One 
needs a large sample of local massive clusters (say a few thousand) to distinguish 
these two models definitely. This conclusion is in good agreement with an early 
study of Jing et al. (1995) on this subject. 
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4. Conclusions 

We have presented a triaxial modeling of the dark matter halo density profiles 
on the basis of the combined analysis of the high-resolution halo simulations (12 
halos with N ~ 10^ particles within their virial radius) and the large cosmologi- 
cal simulations (5 realizations with N = 512^ particles in a 100/i~^Mpc boxsize). 
In particular, we found that the universal density profile discovered by NFW in 
the spherical model can be also generalized to our triaxial model description. 
Our triaxial density profile is specified by the concentration parameter Cg and 
the scaling radius Rq (or the virial radius Rg in the triaxial modeling) as well as 
the axis ratios a/c and a/b. 

We have obtained accurate fitting formulae for those parameters which are 
of practical importance in exploring the theoretical and observational conse- 
quences of our triaxial model. Because the page limit of the proceedings, we 
could not have typed in these formulae. We refer the interested reader to the 
journal paper of Jing & Suto (2002) for the fitting formulae. As Springel sum- 
marized in his talk (this volume), the shape distribution of dark halos is now 
well determined in CDM models. 

The model predictions are compared with various observations of halo 
shapes around galaxies or clusters of galaxies. While a good agreement is found 
between the concordance model prediction and the available observations, a large 
sample of galactic and cluster halos is needed to definitely test the prediction of 
the concordance model or to distinguish among theories of different dark mat- 
ter, because the shape distribution of halos are generally very broad. It is also 
important to point out that iso-potential surfaces are much rounder than the 
iso-density surfaces, especially at halo outer part r ^ Vmr. 
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